Perspectives on Intracluster Enrichment and the Stellar Initial 
Mass Function in Elliptical Galaxies 

ABSTRACT 

The amount of metals in the Intracluster Medium (ICM) in rich galaxy clus- 
ters exceeds that expected based on the observed stellar population by a large fac- 
tor. We quantify this discrepancy - which we term the “cluster elemental abun- 
dance paradox” - and investigate the required properties of the ICM-enriching 
population. The necessary enhancement in metal enrichment may, in principle, 
originate in the observed stellar population if a larger fraction of stars in the 
supernova-progenitor mass range form from an initial mass function (IMF) that 
is either bottom-light or top-heavy, with the latter in some conflict with observed 
ICM abundance ratios. Other alternatives that imply more modest revisions to 
the IMF, mass return and remnant fractions, and primordial fraction, posit an 
increase in the fraction of 3 — 8 M 0 stars that explode as SNIa or assume that 
there are more stars than conventionally thought - although the latter implies 
a high star formation efficiency. We discuss the feasibility of these various so- 
lutions and the implications for the diversity of star formation, the process of 
elliptical galaxy formation, and the nature of this “hidden” source of ICM metal 
enrichment in light of recent evidence of an elliptical galaxy IMF that, because 
it is skewed to low masses, deepens the paradox. 


1. Context 

The hot plasma that pervades the volume of the most massive galaxy clusters - the 
intracluster medium (ICM) - provides a wealth of diagnostic data on the process of galaxy 
formation in structures formed from the largest primordial density fluctuations to have en- 
tered the nonlinear regime and undergone gravitational collapse. The history and efficiency 
of star formation, and the effects of interactions among galaxies and between galaxies and 
the environment in the form of infall, outflow, and dynamical stripping, are reflected in the 
thermal and chemical properties of the ICM - both in individual systems and in the evolving 
population of clusters. 

In massive (> 10 14 M 0 ) clusters, the ICM dominates the baryon inventory and accounts 
for > 10% of the total matter content (e.g., Lagana et al. 2011). The discovery that cluster 
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gas fractions expected to be representative of the universe as a whole exceed flbaryon 1 precip- 
itated a crisis referred to as the “baryon catastrophe” when combined with the assumption 
of a matter- dominated universe and evidence for, and the inflationary prediction of, a flat 
universe - i.e. fl ma tter = 1 (Fabian 1991; White et al. 1993). Of course this paradox was 
resolved by the discovery of dark energy and the concordance cosmology which reconciles a 
flat universe with a reduced D matter . The cosmic baryon matter fraction is now accurately 
determined, D baryon /f2 matter = 0.17 (Jarosik et al. 2011), with a relatively modest percentage 
of the baryons collapsed into stars, fl stars /^matter ~ 0.0074 — 0.011 (Gallazzi et al. 2008). 

ffowever a related paradox persists to this day. While a solar nucleosynthetic yield 
(i.e., star formation ultimately producing the amount of metals necessary to enrich one solar 
mass to solar abundances) is sufficient only to enrich baryons to average abundance of < 0.1 
on a universal scale, cluster baryons are enriched in Fe (and other elements) to ~half-solar 
(Tarnura et al. 2004; de Plaa et al. 2007; Leccardi & Molendi 2008; Bregman et al. 2010; 
Matsushita 2011; Baldi et al. 2012; Andreon 2012). 2 Star formation in cluster galaxies is 
evidently more efficient than in the held; however, a large discordance - first emphasized 
following the groundbreaking accuracy and range of abundance measurements made with 
the ASCA X-ray Observatory - remains between cluster metals and the number of stars 
evidently available to produce these metals (Loewenstein & Mushotzky 1996; Mushotzky & 
Loewenstein 1997; Pagel 1999, 2002; Lin et al. 2003; Finoguenov et al. 2003; Portinari et al. 
2004; Lin & Mohr 2004; De Lucia et al. 2004; Loewenstein 2006; Maoz et al. 2010; Bregman 
et al. 2010). This may be framed in terms of the Fe-mass-to- light ratio (Arnaud et al. 1992; 
Renzini et al. 1993): for M*/ Lb = 5 a solar yield corresponds to M^e/L# ~ 0.0065 - falling 
short by a factor of 5 or more compared to what is measured (Sakuma et al. 2011; Sato et 
al. 2012). 

To resolve this “cluster elemental abundance paradox” one generally must conclude that 
there was more star formation in clusters than conventionally estimated, and/or that star 
formation in galaxy clusters has an enhanced efficiency of producing supernova progenitors 
and synthesizing metals. “Missing” stars may take the form of low surface brightness intra- 
cluster light (ICL), inferred both from observations (Zaritsky et al. 2004; Lin & Mohr 2004; 
Gonzalez et al. 2007) and simulations (Puchwcine al. 2010; Rudick et al. 2011). Recent 
literature includes a large range in calculated star-to-ICM ratios primarily due to divergent 
ICL estimates, but also to different assumed stellar mass-to-light ratios. A stellar initial 
mass function (IMF) that is relatively top heavy increases both the ratio of stars formed, 


1 P represents mean densities in units of the critical density that delineates closed and open universes. 

2 We adopt the solar abundance standard of Asplund et al. (2009) where Fe/H= 3.16 x 10 -5 by number. 
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and of metals created, to present-day light. Such an IMF may be bimodal in nature (Elbaz 
et al. 1995; Larson 1998; Moretti et al. 2003), the second mode perhaps associated with a 
distinct pre-enrichment population (Bregman et ah 2010) where Population III hypernovae 
(Loewenstein 2001) may play a role. Given the conventional wisdom that most ICM metals 
originate in elliptical galaxies (Arnaud et ah 1992), this problem clearly connects to the fun- 
damental galaxy formation questions of the IMF in ellipticals and the transport of material 
from these galaxies into intergalactic space. 

In this paper we undertake a fresh and comprehensive, though generic, examination 
of the metal inventory in the ICM of rich galaxy clusters. We focus on addressing a single 
well-defined, though multifaceted, question: what characteristics of the stellar population are 
necessary to produce the observed level of ICM enrichment? In doing so we address issues 
related to the IMF, star and galaxy formation efficiency, galactic winds, the astrophysics of 
supernova progenitors and explosions, and the apportionment of products of different super- 
nova types into stars and ICM. Section 2 quantitatively summarizes the cluster elemental 
abundance paradox. 

In Section 3, where a standard IMF is assumed, the level of enrichment and abundance 
pattern are related to phenomenological parameters that encapsulate the star formation 
efficiency, and the demographics of supernovae and the success of stars in locking up the 
products of their explosion. Substantial departures from standard values are required to 
match observations. In Section 4, we cast a wider net by considering the effects on ICM 
enrichment of a wide range of IMFs in the context of a self-consistent galaxy chemical 
evolution treatment that accounts for the relevant astrophysical quantities. When juxtaposed 
with recent evidence for an IMF in elliptical galaxies that produces fewer metals than a local 
IMF, this analysis reinforces and clarifies the conflict between the ICM metallicity and the 
characteristics of the stellar population generally assumed responsible for ICM enrichment. 
Results and their implications are discussed, and conclusions summarized, in Section 5. 


2. The Cluster Elemental Abundance Paradox Quantified 

2.1. Basics 

Consider the total baryon mass in a cluster of galaxies within some sufficiently large 
radius that it may be considered a closed box in the chemical evolution sense - that is, all 
products resulting from the transformation of some of this gas into stars (including the stars 
themselves) are contained within this radius. There is direct evidence that this is a good 
approximation for sufficiently massive clusters if invoked at radii that are a significant fraction 



4 


of the virial radius, based on the consistency between the total cluster baryon fraction and 
the universal value mentioned above (Landry et ah 2012), perhaps with a 10% “depletion 1 ' 
correction at r 50 o (the radius within which the average mass density is 500 times the critical 
density); see Gonzalez et ah (2007); Pratt et ah (2009); Giodini et al. (2009); Ade et al. 
(2012); Planellcs et ah (2012). Presumably this is a result of the extreme depth of their 
gravitational potential wells. However it is possible that this also applies to galaxy groups, 
and perhaps even giant elliptical galaxies if one could inventory the gas all the way out to the 
virial radius (and perhaps beyond, if entropy injection has dispersed the gas distribution). 

We define the overall efficiency of converting gas into stars, £ s f, such that the total mass 
in stars formed (regardless of where) is 




*,/orm 


= £ o f Ml 


sf ■LVJ-baryom 


(i) 


where Mb aryon is the total baryon mass being considered. At the present time the total mass 
in stars, whether contained in individual cluster galaxies (including the brightest cluster 
galaxy - BCG) or associated with intracluster light (ICL), is 

A/* M^baryon ^ sf (f ^*)> (2) 


and the mass in gas is 


M gas 


Mbaryon M* 

Mbaryon [1 £s/(l G)] > 


( 3 ) 

( 4 ) 


where r* is the stellar “mass return fraction” - the fraction of the mass previously formed into 
stars recycled back into gas. For massive clusters one may neglect the distinction between 
the total mass in gas and the mass in the ICM, and henceforth we equate M gas with Mjcm- 
The star formation efficiency, in terms of the observable Micm/M * is 


£ sf 


(1 — r*) 1 



Mjcm \ 1 
M* ) 


( 5 ) 


We consider the enrichment of cluster baryons in chemical elements released in supernova 
explosions, i.e. those of atomic number A > 8. It is straightforward to extend the analysis 
to elements statically synthesized in intermediate mass stars; however, the elements (C, N) 
to which this applies are not well-constrained by current X-ray observations. 

Both the overall level of baryon enrichment (that is, the metallicity) and the abundance 
pattern are determined by the total number of supernovae and their nncleosynthetic yields. 
We separate the enrichment contributions of the two main classes of supernovae - Type 
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la supernovae (SNIa) that result from the explosion of a white dwarf, and core collapse 
supernovae (SNcc). Their total numbers may be expressed as 


N cc 


V M * t form 


( 6 ) 


and 

N Ia = V Ia M^ form , (7) 

where rf c and r] Ia are, respectively, the specific numbers of SNcc and SNIa explosions per 
star formed. It is useful to define the total supernova number, supernova ratio, and SNIa 
fraction as follows: 



N sn = N cc + N Ia 

(8) 


= V SN (1 - 

(9) 

where 

V SN = ^cc + 7] l a . 

(10) 


r sn = 

rf c ’ 

(11) 

and 

la tdSN 

cl a _ 11 

1 ~ V SN 1 + R SN ' 

(12) 

Note that the number of SN per unit mass in the ICM is 



" SN SN 1 1 ,-l /• 

M, CM =ri (1 “ r>) W 

(13) 

where /* and f IC M 

are the present-day mass fractions of stars and gas: /+ / frc'M 

= M*/M IC m ■ 


2.2. Stars and Supernovae 

A combination of theoretical and empirical considerations enter into determination of 
the stellar and supernovae parameters. Star formation is not sufficiently well-understood to 
allow an a priori estimate of e s f, and M* is estimated from observations of cluster starlight 
in galaxies and intracluster space. To convert to mass, stellar population synthesis can be 
employed, involving assumptions about the IMF - among many other factors. The rnass- 
to-light ratio in individual galaxies can also be inferred from dynamical modeling of stellar 
velocity dispersion distributions. The mass return fraction, r*, may be calculated from the 
IMF, star formation history (SFH), and relation of stellar remnant mass to progenitor mass 
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for individual stars estimated from stellar evolution theory and observations in the Galaxy. 
The SNcc efficiency, rj cc , depends on the IMF and range of masses that result in core collapse 
explosions. Although one may model the time-dependence of the SNIa rate from assumptions 
about the binary star progenitor configuration and the distributions of binary mass ratios 
and separations, the normalization is difficult to estimate a priori. Empirical estimates of y Ia 
(and rf c ) are emerging from supernova surveys that also constrain the distribution of delay 
times from the observed evolution of the SNIa rate (Maoz & Mannucci 2012; Sand et al. 
2012). Estimation of these fundamental quantities - M*, r*, rj Ia , rf c - in principle requires a 
convolution and a time-integration over an evolving galaxy population with disparate SFHs 
and, conceivably, IMFs. A first order approximation treats the ensemble stellar population as 
a single simple population of stars formed at some (high) redshift with a common IMF - an 
approximation most suitable for rich clusters where the total galaxy mass is most dominated 
by elliptical galaxies. 

We can insert some reasonable values to get a sense for the expected level of super- 
nova enrichment and relative contribution from the two classes of supernovae. For a “diet 
Salpeter” IMF that represents a simple alteration - proposed as a means of reconciliation 
with the observed relative frequency of ~subsolar mass stars (Bell & de Jong 2001) - of 
the classic single-slope Salpeter function (Salpeter 1955), rj Ia ~ 0.002 (Maoz & Mannucci 
2012), r* ~ 0.35 for an old stellar population (Fardal et al. 2007; O’Rourke et al. 2011), and 
y cc ~ 0.008 (Maoz & Mannucci 2012; Botticella et al. 2012; Dahlen et al. 2012). 

One then predicts R SN ~ 0. 25 (/*• rsj 0.2) and r/ SN ~ 0.01. For massive clusters (i.e., 
M 5 oo = M(r< 500 ) > 10 14 M 0 ), recent studies report a range in stellar mass fraction evaluated 
at r 5 oo, reflecting different treatments of ICL and in conversion from light to mass (Zhang et 
al. 2011), with /*/ ficM typically ~ 0.1 (Lin et al. 2003; Gonzalez et al. 2007; Giodini et al. 
2009; Ettori et al. 2009; Dai et al. 2010; Andreon 2010; Bregman et al. 2010; Lagana et al. 
2011; Balogh et al. 2011; Zhang et al. 2011; Lin et al. 2012) (though with substantial system- 
atic uncertainty and study-to-study variation; see Leauthaud et al. (2012)), and evidence 
of an increase in magnitude and scatter with decreasing cluster mass. The resulting total 
number of supernova explosions per solar mass of ICM is ~ 1.5 x 10 3 (10 f*/ J icm) Mq 1 . 


2.3. Application to the ICM 

The equation for the mass of the ith element in the ICM, Mj , in terms of the number 
of SNIa and SNcc that enrich the ICM, N Ia and N cc , and the yields per SNIa and SNcc, y\ a 
and (yf ) , is 


Mi = N cc (y f) + N Ia yj a = N sn ( 1 + A SiV )- 1 ((yf) + R SN yl a ), 


(14) 
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where, as defined above, R SN = N Ia /N cc and N SN = N Ia + N cc . The IMF ( d)-averaged 
SNcc yield is 

f' nup dm<?>(m)yf c (m) 

I cc\ _ Jmcc V > /i 

\i>i / rniup 7 i/ \ i \- L ' J ) 

J m amcpyrn) 


where rn cc and rn up are the lower and upper limits for the masses of SNcc progenitors, and a 
single universal set of SNIa yields is assumed. Despite a plethora of IMF parameterizations 
(Kroupa et al. 2012), there is general consensus that a Salpeter slope (Salpeter 1955), f ~ 
m ~ 2 ' 35 applies at the high mass end relevant for SNcc - at least for star formation under 
“normal” conditions. 


The resulting mass fraction in the ICM (mass Mjcm ) of the ith element, /, , is 


M. n sn 

h = (i + R SN )-\{y? ) + R SN y! a ), 

m icm m tcm 

and the mass fraction relative to the solar mass fraction, 

fj _ N SN ( 1 , E»SJV\-1 f fyT) , D SN yl a 

fi 0 M icm [ + ’ \ fie he 


(16) 


(17) 


The relationship between mass fraction, fi, and abundance, z t , of the ith element (the number 
of atoms of element i relative to that of H, i.e. the entries in standard abundance tables) is 
z % = (fa/ X)(Ai/ A h)- 1 , where X and A H are the hydrogen mass fraction and atomic weight 
(Ah = 1.008 AMU) and A i: the atomic weight of the ith element. Relative to solar, the 
abundance is 

a = A 

z ie X he fie 

where /j Q = Zi & X & (Ai/ Ah), and the approximation X = X & is invoked - an approximation 
that is valid as long as the total mass fraction of metals (< 2% for solar abundances) is 
small and the He abundance is fixed. X for various solar standard abundance sets is given 
in Table 4 of Asplund et al. (2009). 



Finally, the abundance relative to solar is expressed as 

,. j\r SN 

Zi = ^- = — (l + R™)- 1 ^?) + R SN y! a ') (19) 

z ie m icm 

where N SN / Mjqm is given by equation (13); and, y = (yf c )/he an< l y\ a ' = A/ he- F° r 
Micm measured in M 0 these are the yields of the ith element relative to the mass of that 
element contained in one M 0 of solar abundance material. As noted above for the total 
baryons, R SN and N SN fully determine the level and pattern of ICM enrichment - modulo 
sets of SNIa and SNcc yields - and can be compared to ICM abundances. The new approach 
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of Bulbul et al. (2012) directly fits X-ray spectra to a model parameterized by R SN and N SN 
via the abundance predicted by equation (19). 

Focusing on Fe, the element with the most widely determined and most accurate global 
ICM abundance measurement, equation (19) predicts Z Fe j C M = 0.255(10 /*///cm)> f° r 
the values of R SN and N sn derived at the end of the previous subsection and adopting 
(yj a , (|/j CC )) = (0.743 M 0 , 0.0825 M 0 ) from Kobayashi et al. (2006) - about half the typical 
observed value for /*/ ficM = 0.1. However this assumes that all of the metals produced by 
supernovae reside in the ICM. The values R SN and N SN relevant here correspond to those 
supernova explosions that enrich the ICM (or, for some particular X-ray measurement, those 
in a particular spectral extraction region of a particular cluster). Not all of the products 
resulting from supernova nucleosynthesis are available to enrich the ICM. 


2.4. Metals Locked Up in Stars 

One approach to evaluating galaxy cluster enrichment is to estimate the total inventory 
of metals in stars and in the ICM in the context of the total required number of supernova 
explosions. However, our focus here will be on the ICM which offers more accurate abundance 
determinations over a wider range of elements via X-ray spectroscopy. This specifically 
requires a correction accounting for how supernova products are apportioned among gas and 
stars. 

The galactic mass in rich clusters is dominated by early-type systems that form their 
stars rapidly. This results in the well-established enhancement in [a/Fe], the abundance 
ratio of a-elements to Fe (expressed as the logarithm with respect to solar) - i.e., SNcc 
products are preferentially locked up in stars. In their investigation of the giant elliptical 
galaxy NGC 4472, Loewenstein & Davis (2010) found that a ratio of SNIa to total supernovae 
of N^/N™ ~ 0.11 (N?/N? n r 0.89) and number of supernova per mass in (present-day) 
stars of N* n /M* ~ 0.0083 resulted in [a/Fe]* ~ 0.25 and Z Fe ~ 1 (as observed in this 
particular galaxy, but typical of the class; see also Lin et al. 2003, Gallazzi et al. 2008). This 
enables us to estimate the lock-up corrections, rff and rjl a , needed to convert the specific 
supernova numbers to those available for enrichment of the ICM as follows: 

N Ia 

Vl a = -^-(i - r *) « 6.0 x 10 ~ A Z Fe ^ (20) 

N cc 

V? = jj-0- ~ u) « 4.8 x 10-W 


and 


( 21 ) 
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The corresponding values available to enrich the ICM, t)jq M and rffcMi are 

vicM = n Ia - vl a = 1.4 X 10- 3 M- 1 , (22) 

and 

vfcu = »“ - V? = 3.2 x 1(T 3 M- 1 (23) 

for the default parameters considered above. That is, the metal production from ~ 60% of 
SNcc and ~ 30% of SNIa must be locked up in stars to enrich them to solar Fe abundances and 
[a/Fe]* ~ 0.25 - with these factors deducted from the total to obtain the effective enrichment 
of the ICM. A relative overabundance of SNIa contributing to ICM enrichment is expected 
based on the inference that the accelerated formation of stars in clusters preferentially locks up 
the products of SNcc. This asymmetry is observed in the abundance patterns in cluster cores 
(de Plaa et al. 2007; Lovisari et al. 2011), although whether this extends globally remains 
an open question - e.g., the smothering of galactic winds in central dominant galaxies may 
skew the pattern, relative to the ICM as a whole, via concentrated direct injection of SNIa. 

Based on these estimates, for the supernovae remaining available to enrich the ICM, 
R sn ~ 0.44 ( f Ia ~ 0.30) and N sn /M icm ~ 0.71(10 f*/ f ICM ) x 10~ 3 Mq" 1 . This enriches 
the ICM in Fe to the level Z Fe ,iCM = 0 . 155 ( 10 /*// jcm ), thus quantifying the paradox that 
baryons in clusters of galaxies are enriched beyond what is expected based on the stars we 
see in galaxies today - unless either the star formation efficiency exceeds that in the field by 
a factor of ~ 3, or supernovae are more efficiently produced per unit star formation. 


3. A More Comprehensive Examination (I) 

In this, and subsequent, sections we investigate stellar and ICM abundance predictions 
for a range of elements - focusing on a subset selected on the basis of a combination of 
accessibility and diagnostic power: O, Mg, Si, Fe, and Ni. First, we consider a wide range 
of published yield sets and apportionment of metals into stars and ICM in an effort to place 
robust constraints on the required efficiency of star formation. In order to be as general 
and assumption- free as possible we introduce two parameters that gauge the efficiency with 
which stars may lock up supernova products and do so asymmetrically, i.e. preferentially 
for SNcc relative to SNIa. We consider the specific effects of varying the IMF, with its 
coupled impact on r*, rf c , and rj Ia , in the following section, fixing these parameters at the 
values described above (r* ~ 0. 35, rj Ia rs_/ 0.002, rj cc ~ 0.008) for immediate purposes. Note 
that with these parameters, supernovae enrich cluster baryons in Fe to a mass-averaged 
abundance of Z barFe = 1.66 e s f. 
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We generalize our treatment of quantifying the fraction of metals synthesized by super- 
novae that are inaccessible for ICM enrichment due to lock-up in stars by defining a SNcc 
lock-up fraction, /3 CC , and supernova asymmetry parameter, a SN , such that 


f3 cc = 


Vi 

Tj c 


(24) 


and 

/ n SN \ -1 

^ (§») ■ < 25 > 

where R* N = r][ a /rff. By definition, /3 CC < 1, and a SN > 1 is expected under conditions 
where rapid conversion of gas into stars results in preferential incorporation into stars of 
SNcc products with respect to those from SNIa products that are released over a relatively 
extended time interval. We provisionally adopt the values that correspond to the estimates 
of the previous section as standard for the remainder of this section: /3 CC = 0.6 and a SN = 
2. ICM abundances are calculated, for a given star formation efficiency e s f, from ICM- 
specihc versions of equations (5), (10), (11), (13), (19), with the supernovae per star formed 
effectively reduced to 

Vicm = V Ia ^1 - -^sN^j , (26) 

and 

Vi c cm = V cc ( l-/n- (27) 

The results, assuming Kobayashi et al. (2006) supernova nucleosynthetic yields, 3 are dis- 
played in Figure 1. Figures la and lb confirm that the default (IMF and) lock-up parameters 
predict a stellar population enriched to solar Fe abundances, and ICM Fe abundances lower 
than observed for e s f ~ 0.15 (/*/ ficM = 0.11). One may recover the observed level of ICM 

Fe enrichment for e s j ~ 0.3 (/*/ ficM = 0.24), but only for extreme values of the lock-up 

parameters that imply that most of the Fe produced by stars ends up in the ICM, and that 
stellar Fe abundances are well below solar - in contradiction to observations. 


Figure lc plots Zp e ,icM and Zp e ,* versus e s f for selected pairs (/ 3 CC , a SN ) - including 
the default. Also shown is the limiting case (3 CC — > 0 where 100% of metals produced by stars 
reside in the ICM. In this case 


ZFe,ICM 


lMe s f 

1 - e sf ( 1 - r*)’ 


(28) 


3 Although these are averages for an IMF with a Salpeter slope at the high mass end, we adopt these in 
general - the differences for the IMFs we consider are generally small. 
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Fig. 1. — Left panel (a): ICM (solid lines) and stellar (broken line) Fe abundance versus lock-up fraction 
/ 3 CC for (the adopted standard) a SN = 2, and e s / = 0.15 or 0.30. The horizontal dotted line shows the 
typical observed value abundance, Zp e ,iCM = 0.5, the vertical dotted line the adopted standard (3 CC = 0.6. 
Middle panel (b): Same as (a) versus a SN for (the adopted standard) (3 CC = 0.6. The vertical dotted line 
shows the adopted standard a SN = 2. Right panel (c): ICM (curves) and stellar (horizontal lines) Fe 
abundance versus e s / for pairs (a SN , /3 cc )=(10,0.1)(short-dashed), (3,0.2)(dotted), (2,0.6)(standard: solid), 
(2. 0,0. 3) (long-dashed), and (l,0.8)(dot-short-dashed). The dot-long-dashed line shows the limiting case 
(3 CC -» 0 . 

from which on can see that e s f > 0.25 represents an absolute lower limit to the star formation 
efficiency required to enrich the ICM to Z Fc j C m > 0.5. This figure provides an alternative 
demonstration that, for e s f = 0.15, Z FcI cm < 0.3 even for extreme models where such a 
large fraction of supernova-produced metals is released into the ICM that insufficient metals 
remain available to enrich the stars to the observed level. Both relatively large star formation, 
and small lock-up efficiencies 4 are required to simultaneously enrich the stars and ICM to 
the observed level (see, also, Sivanandam et al. 2009). 

The increasing divergence of stellar and ICM abundance ratios (that are independent of 
e s f) with increasing a SN is shown for /3 CC = 0.6 in Figure 2a. One can see how, in this case, 
the enhanced [a/Fe] measured in the old stellar populations that dominate cluster galaxies 
implies a large asymmetry parameter and, as a result, subsolar ratios of a-elements with 
respect to Fe in the ICM. Since [ a/ Fe\icM << 1 is not observed, large values of a SN may be 
ruled out. This divergence narrows with decreasing f3 cc and ~solar ICM abundance ratios 
emerge for /3 CC ~ 0.3, i.e a lower lock-up fraction (Figure 2b). 

The effects on cluster enrichment of adopting different SNIa yield sets is shown in Figures 
3-5. The solid lines correspond to previous plots that utilize Kobayashi et al. (2006) yields, 
the dotted lines use the same SNcc yields, but alternative SNIa yield sets from Nomoto et 


4 or, more precisely for Fe, (3 cc /a SN « 1 
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Fig. 2. — Left panel (a): ICM (solid lines) and stellar (broken lines) abundance ratios with respect to Fe 
for O (red), Mg (blue), Si (green), and Ni (yellow) versus a SN for the standard f3 cc = 0.6. Right panel (b): 
same as (a) for /3 CC = 0.3. 


al. (1997); Maeda et al. (2010) - see Loewenstcin & Davis (2010). Here we focus on Fe, Ni, 
and Si (O and Mg are always dominated by SNcc enrichment and insensitive to this choice). 
Figure 3 shows the total baryon enrichment, which is independent of f3 cc and a SN , and 
demonstrates that star formation efficiencies £ s f ~ 0.25 — 0.5 are required to enrich cluster 
baryons to a relatively modest Fe abundance of half-solar - e s f ~ 0.4 — 0.7 to attain 0.75 
solar. Figures 4 and 5 that, respectively, show the ICM abundances for the standard lock-up 
parameters, and the maximum ICM abundances corresponding to zero stellar metallicity, 
demonstrate that the conclusions about the required efficiency of star formation are not a 
result of a particular choice of yields sets. Using SNcc yields from Woosley & Weaver (1995) 
(their standard explosion energy, solar abundance progenitor model) does not alter these 
conclusions. 

From the results in this section, we confirm that ICM Fe abundances cannot be produced 
if £ s f ~ 0.15 (/*/ ficM ~ 0.11) and quantified the shortfall as a function of how efficiently 
metals in general, and SNIa products in particular, are locked up in stars. For £ s f ~ 0.3 
(/*/ f ICM rs./ 0.24), they can - but only for small lock-up fractions such that ~ 85% of 
SNIa, and ~ 70% of SNcc, metal production is embedded in the ICM and (by implication) 
%Fe,* ~ 0-5- The assumption of a standard IMF is adopted throughout this section, an 
assumption we relax in the following section. 
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Fig. 3. — Left panel (a): Average Fe abundance of all cluster baryons (stars and gas) as a function of star 
formation efficiency, e a f, assuming a single SNcc yield set but a wide range of SNIa yield sets (see text for 
details). The horizontal dotted line shows the very conservative value of Zp e ,b or = 0.5. Middle panel (b): 
Same as (a) for Ni. Right panel (c): Same as (a) for Si. 
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Fig. 4. — Same as Figure 3 for ICM abundance, assuming the standard (3 CC = 0.6, a SN = 2. 


4. A More Comprehensive Examination (II): Effects of Changing the IMF 

The initial mass function (IMF) of stars in cluster galaxies and intracluster space is 
intimately connected to estimates of ICM enrichment, impacting the mass in stars calculated 
from the total light, the mass return fraction (or, equivalently, the ratio of current stellar 
mass to mass converted into stars), and the numbers of SNcc and SNIa explosions expected 
per mass formed into stars. The general characteristics of the IMF in various Milky Way sub- 
populations are now well-determined and generally mutually consistent (Bastian et al. 2010). 
Although several functional forms are commonly used for the “canonical” IMF (Kroupa et 
al. 2012), these must share the properties of a Salpeter-like slope at high mass with a break 
to a flatter slope below ~ 0.5 — 1 M 0 . Many of the best-studied environments are consistent 
with the hypothesis that this form is “universal” in space and time, but variations in extreme 
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Fig. 5. — Same as Figure 3 for the maximum ICM abundance, i.e. (3 CC = 0. 

star forming environments that predominate in the early universe (when most stars in the 
elliptical galaxies that dominate the stellar content in clusters form) may explain a number of 
anomalies, such as an apparent inconsistency between the observed evolution of the global 
star formation rate and stellar mass densities (Narayanan & Dave 2012a, and references 
therein). 

Several recent observational investigations of elliptical galaxies find direct evidence in 
the mass-to-light ratio for either an excess of stellar remnants as realized in a “top-heavy” 
IMF; or, of low mass stars as realized for a “bottom-heavy” IMF (e.g., Cappellari et al. 
2012a; see Section 5 below). By exploiting the level and pattern of ICM abundances we 
constrain the properties of the enriching stellar population. We may then exclude particular 
elliptical galaxy IMFs under the parsimonious assumption that these optically studied stars 
originate from the same parent IMF as those that enrich the ICM - or, alternatively, call 
this assumption into question. 

4.1. Models and Parameters 

Simply put, the level of metal enrichment of the stellar and ICM baryonic sub-components 
in clusters is a reflection of their respective total masses, the total numbers of SNIa and SNcc 
that enrich each constituent, and the nucleosynthetic yields of each of these supernova explo- 
sion types. In previous sections, these are expressed in terms of the mass return fraction, r* , 
and formation efficiency, e s f of the stars, the total and relative numbers of SNIa and SNcc 
per star formed ( R SN , f Ia ), and phenomenological supernova lock-up and asymmetry pa- 
rameters (/3 CC and a SN ) characterizing the ultimate destination (ICM or stars) of supernova 
products. In Appendix A these are further deconstructed into more fundamental astro- 
physical functions and parameters directly connected to stellar and galaxy evolution, thus 
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providing a self-consistent astrophysics! framework for understanding the effect of varying 
the IMF on cluster enrichment. Ultimately these are reduced to the following: (1) The func- 
tional form of the IMF (see below); (2) the present-day main sequence turnoff mass (0.9 M 0 ); 

(3) the remnant-progenitor mass relationships for (< 8 M 0 ) intermediate- (equation A4) and 
(> 8 M 0 ) high-mass (equation A29) stars derived, for the former, from well-established stan- 
dard white dwarf masses and, for the latter, from a model for the evolution of the stellar 
population and the delayed-explosion compact remnant prescription in Fryer et al. (2012); 

(4) the ratio of mass ejected from galaxies into the ICM during star formation to the mass 
of stars formed, 5gw] (5) the galaxy formation efficiency, e ga i ; 5 and, (5) various supernovae 
switches and parameters that we now describe. For SNcc we assume progenitor masses from 
m cc = 8 M 0 to m up , where m up may differ from the IMF upper mass limit rrihi (but is the 
same by default, and assumed so in calculating the high mass return fraction). Since we 
adopt the IMF-averaged SNcc yields of Kobayashi et al. (2006) as a function of progenitor 
metallicity Z cc , Z cc must be specified as well. For SNIa we consider the yield sets described 
in Section 3, assume progenitors in the 3 — 8 M 0 range, and must specify the efficiency e Ia 
defined as the fraction of 3 — 8 M 0 that result in SNIa. In addition, the “prompt” fraction 
of SNIa that explode during the star formation epoch and so may be incorporated into stars 
or ICM, /p“, (while a fraction 1 — / strictly enrich the ICM) must be specified. 

Our default IMF is the Kroupa et al. (2012) segmented power-law with slopes and mass 
scales that can explain local star formation (Appendix A); other defaults are Z cc = 1, W7 
yields, e ga i = 0.25 and dew — 0.5 (e s f = 0.17, /*/ ficM ~ 0.11), e Ia = 0.076 (Section A. 2), 
and /p“ = 0.5. Under these conditions, r* = 0.41 (r* m = 0.25, = 0.16 for intermediate- 

and high- mass stars, as delineated above) while the fraction in stellar remnants is 0.17 
(0.11/0.06 from intermediate/high- mass stars), r] cc = 0.011, rj Ia = 0.0022, a SN = 2 and 
(3 CC = 0.39 - similar to the default parameters in Section 3. The resulting abundances are 
shown in Table 1 where, once again as expected, we find reasonable stellar abundances but 
ICM abundances too low by a factor of ~ 2. For comparison we also display the abundances 
for a Salpeter IMF, which fails to provide sufficient metals for all components (including 
stars) for this default set of parameters. 6 

Our approach to examining the effects of varying the IMF on ICM enrichment, that 
attempts to make comparisons at fixed values of the observables to the extent possible, is 
as follows. As detailed in Appendix A, we consider departures from the “canonical” IMF 


5 Defined as the fraction of baryons initially in galaxies, this is related to the star formation efficiency 
defined in equation (2) by the expression e s f = £ so z(l + 5gw)~ 1 ■ 

6 It should be noted that most observational estimates of /*///CM are IMF-dependent, and would be 
larger for a Salpeter IMF given a fixed amount of optical light. 



Table 1. ICM Abundances from Default Parameters 



0 

Mg 

Si 

Fe 

Ni 

Baryons 

Stars 

ICM 

0.44 (0.27) 
1.74 (1.05) 
0.30 (0.15) 

0.30 (0.18) 
1.19 (0.72) 
0.20 (0.10) 

0.38 (0.23) 
1.35 (0.82) 
0.27 (0.15) 

0.32 (0.20) 
0.87 (0.54) 
0.26 (0.15) 

0.84 (0.53) 
1.90 (1.21) 
0.72 (0.43) 


Note. — All abundances relative to (Asplund et al. 2009) solar standard. The 
values in parentheses are for a Salpeter IMF with the default values of e ga i , Sqw ■> 
e Ia , and fp a , and the same range of masses (0.07 — 150 Mq). 


Table 2. Model IMF Parameters 
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0.07 


1.8 -»■ 2.3 
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sl-1 

0.07 

0.5 

0.3 2.4 

2.3 

sl-2 

0.07 

1.0 

0 
CO 

1 

1 ° 

CO 

2.3 

sl-3 

0.07 

0.5 

1.3 

1.8 -> 2.5 

sl-4 

0.07 

1.0 

1.3 

10 

CN 

t 

00 

si- 5 

0.07 

8.0 

1.3 

1 

w 

CO 

m-1 

0.01 -> 1.0 

1.0 

1.3 

2.3 

m-2 

0.01 -t- 0.5 

0.5 

1.3 

2.3 

m-3 

0.07 

0.5 -*• 5.45 

1.3 

2.3 


Note. — irihi = 150 Mq for all models displayed here. 
m 6r — equivalent to m2 = 0.5 M© for models sl-1, sl-3 , and 
m-2 ; and m3 (default: 1 M©) for models sl-2 , sl~4, si- 5, m-1 , 
and m-3 — is defined as the mass where the IMF slope transi- 
tions to its high-mass value. Models si- 4 and si- 5 include an 
additional break from ai =1.3 to 0:2 = 2.3 at m2. ot\ = ai 0 
and 0:3 = OLhi, while 0:2 = oli 0 for models sl-2 , m-1, and 
m-3 ; 0:2 = 0:3 for models sl-1, sl-3, m-2; and, is set at the 
default 0:2 = 2.3 for models sl-1, sl~4, si- 5, and m-2. Param- 
eter ranges correspond to those with physical solutions; i.e., 
fp a and Sqw > 0 and e ga i < 1 - for the sl~4 model with 
higher SNIa efficiency (star-to-gas ratio) the range shifts to 
1.5 — ► 2.5 (2 — ► 2.5); see below and Tables 3 and 4. 
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(Table 2) described according to sequences with either a single slope below the break mass 
m3, or an additional distinct slope below 0.5 M e . Either the lower mass limit (m/ G ), or the 
slopes (o;i or a 3 ) at low or high mass ends may be varied. Additionally, m 3 may vary in 
sequences of the first type. To isolate the effect of the IMF on ICM abundances, we impose 
invariance on the stellar Fe and O abundances (Table 1). This determines the necessary 
adjustments in the parameters Sow and f£ a (for some e Ia ). For fixed yield sets, invariance 
in Zp e ^ and Zq,* assures invariance in all stellar abundances. Finally, we consider these 
variations at fixed present- day baryon inventory (11% stars, 89% ICM), which is equivalent 
to adjusting e ga i so as to maintain constant e s f( 1 — r*) (equations 2, A26) - thus enabling 
us to investigate what adjustments in IMF (if any) may explain observed ICM abundances 
assuming this nominal star-to-gas ratio. The imposition of these constraints rule out those 
IMFs that imply unphysical values of $gw{< 0), //“(> 1), or £ ga i{> 1) - i.e., some IMFs are 
incompatible with observed stellar abundances and a ~9:1 ratio of ICM to stars (thus the 
limited range of the variable IMF parameters in Table 2). 


4.2. Impact of Varying the IMF 

We remind the reader that our standard IMF has slope a = 1.3 below 0.5 M 0 , and 
2.3 above. Figures 6a-d show the impact on ICM enrichment of varying one (and only one) 
of the slopes and (in some cases) adjusting the break mass, by plotting the Fe abundance 
and Mg/Fe, Si/Fe, and Ni/Fe ratios versus the deviation in the non-fixed slope from these 
standard values. 7 This covers many of the IMFs considered in the literature as possibly 
resolving various conflicts between expectations and observations of stellar populations in 
elliptical and/or starburst galaxies. For models sl-1 ( sl-2 ) a below 0.5 M 0 (1 M 0 ) is varied 
with a above these single break masses fixed at 2.3 - i.e. positive (negative) Aa corresponds 
to bottom-heavy (-light) IMFs. For models sl-3, sl-4, and sl-5, a is varied above break 
masses of 0.5 M 0 , 1 M 0 , and 8 M 0 , respectively, with a = 2.3 between 0.5 M 0 and the 
break mass in the latter two. For these models, positive (negative) Aa corresponds to top- 
light (-heavy) IMFs. In addition, we plot the results for a single-slope IMF ( sl-0 ). We can 
see that Zp e jcM ~ 0.5 is predicted for either a (1) bottom-light IMF with a(< 1M 0 ) ~ 1, 
(2) top-heavy IMF with a(> O.5M 0 ) ~ 2 or a(> 1M 0 ) ~ 1.8, (3) single-slope IMF (i.e. both 
bottom- and top-heavy) with a ~ 1.8. As expected, Zp e jCM ~ 0.1 for the pure Salpeter 
IMF (single slope and bottom- light cases with slope 2.35). 


'Mg is almost exclusively synthesized in SNcc, Si primarily (but not exclusively) in SNcc, Fe in both 
SNcc and SNIa, and Ni primarily in SNIa. 
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Varying the lower mass cutoff provides an alternative means of producing either bottom- 
heavy (■ nrii 0 < 0.07 M 0 ) or bottom-light (m; G > 0.07 M 0 ) IMFs, while increasing the break 
mass that delineates a = 1.3 from a = 2.3 (Narayanan & Dave 2012a) is also bottom light 
in the sense that the fraction of low-mass stars relative to those at intermediate and high 
mass is suppressed. As shown in Figure 7a, these alternatives can explain Zp e jCM = 0.5 for 
IMFs with the standard high mass slope a = 2.3 if the IMF lower mass cutoff is shifted from 
the standard = 0.07 M 0 to m/ 0 ~ 0.2 M 0 for an IMF with a single break at 1 M 0 , what 
is (essentially) a Salpeter IMF with mi 0 ~ 0.5 M 0 , or for an IMF with the break between 
a = 1.3 and a = 2.3 shifted from 0.5 to ~ 1 M 0 . The last would seem to represent a 
particularly modest departure from the standard IMF. 

Bottom-light and top-heavy scenarios may be directly distinguished in ICM spectra via 
abundance ratios, as demonstrated in Figures 6b-d, 7b-d, and 8a-c; and Table 3. These 
essentially define two branches in the Zp e ,iCM — (a/ Fe)icM plane (Figure 8), with ratios 
in the top-heavy branch connecting to the stellar ratio as the IMF flattens (and f , — > 1), 
and abundances for the bottom-light IMFs (assured to have f,p a = 0.5) universally rising 
in lockstep with the elimination of low-mass stars. Table 3, confined to those models that 
predict Z Fe jcM = 0.5, illustrates how the a/ Fe Fe ,iCM ratios might be exploited to distin- 
guish among bottom-light and top-heavy IMF explanations for ICM enrichment - for the 
former, ratios deviate more strongly from those in stars (smaller a/ Fe Fe jcM, larger Ni/Fe) 
and provide a better match to the data (Simionescu et al. 2009). 


4.3. Implications of Models with Nonstandard IMFs 

Table 4 displays the essential characteristics of models that produce Z Fc jqm = 0.5, and 
are constrained to match the standard values of stellar metallicity and /*/ ficM - 8 Several 
general properties, as well as others that distinguish top-heavy from bottom-light solutions 
emerge. Relative to the model with canonical IMF, all have a relative deficiency of unevolved 
low-mass stars, and hence ~ 30% higher mass return (r*) and remnant (/ rem ) fractions - 
implying ~ 30% upward adjustments in the integrated mass of stars formed based on the 
present-day mass, and ~ 60% upward adjustments based on the luminous stellar mass. The 
“extra” metals are explained by a larger fraction of stars in the supernova-progenitor mass 
range, and a larger ratio of mass in stars formed to present-day stellar mass. 

Successful top-heavy models are characterized by a large “prompt” fraction of SNIa and 
prodigious galactic winds, and have low lock-up fraction and modest asymmetry between 


8 Effects of relaxing the latter are discussed shortly. 
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stellar and ICM abundance patterns. That is, the extra ICM metals are mostly associated 
with the rapid early star formation epoch where both formation of SNcc and SNIa progeni- 
tors, and delivery of the metals from star forming sites to extragalactic hot gas, are efficiently 
realized. 

Successful bottom-light models are characterized by the standard prompt SNIa fraction 
and more modest (though still substantial) galactic winds, with a larger fraction of ICM 
enrichment occurring during the passive, post-star-formation phase. In these models, the 
implied fraction of mass originally in galaxies is 0.42 and the fraction of the ICM that is 
primordial, f pr im = 0.65 - compared to 0.25 and 0.83 in the canonical model. For the top- 
heavy models these take on more extreme values - each on the order 0.5 - although, due to 
the effects of galactic winds, the star formation efficiencies are not appreciably different. 

The different abundance patterns predicted in bottom-light and top-heavy models is 
reflected in values of Rjcm that are ~twice as high for the former, more consistent with - 
though still lower than - the value recently inferred by (Bulbul et al. 2012) for the central 
region in Abell 3112. 


4.4. Other Variations 

4 . 4 .I. Supernova Yields 

In most cases, varying the SNIa yield set primarily affects the predicted (stellar and 
ICM) Ni/Fe ratios - which are not well-determined in clusters at this time. Exceptions are 
yield sets with particularly low (< 0.4 M 0 ) Fe yields, i.e. the C-DEF and C-DDT in Maeda 
et al. (2010). These models cannot self-consistently produce the observed stellar abundances 
and Z Fc jcm = 0.5. Since there is a narrow range of Fe yields in SNcc calculations, our 
results are insensitive to the choice of Z cc - though, in principle, abundance ratios among 
a-elements could carry SNcc yield diagnostic information. Decreasing the SNcc upper mass 
limit, m up , with respect to the IMF upper limit, m^, lowers the predicted metallicities - 
though the effect is small unless the upper mass IMF slope is very flat. 

Overall our models are conservative in the sense of maximizing Fe yields by adopting 
high values for the SNIa Fe yield (0.74 M 0 ), and for rn hl (150 M 0 ). 


4- 4- SNIa Efficiency 

Maoz & Mannucci (2012) estimated that a 1.7x higher Type la supernova rate per unit 



Table 3. Abundance Ratios in Models with Zp e jCM = 0.5 





O/Fe 

Mg/ Fe 

Si/Fe 

Ni/Fe 

stars 



2.0 

1.4 

1.6 

2.2 

salpeter 



1.0 

0.68 

0.97 

2.8 

canonical 



1.1 

0.77 

1.0 

2.7 

sl-0 

a = 

= 1.83 

1.9 

1.3 

1.5 

2.2 

si- 2 

Oilo 

= 1.06 

1.2 

0.85 

1.1 

2.7 

sl-3 

CX-hi 

= 1.97 

1.7 

1.2 

1.4 

2.4 

{e Ia = 0.13) 

CX-hi 

= 2.22 

0.85 

0.57 

0.87 

2.9 

(star/gas=0.25) 

&hi 

= 2.36 

1.0 

0.69 

0.97 

2.8 

sl-4 

CX-hi 

= 1.85 

1.9 

1.3 

1.5 

2.3 

m-1 

mio 

= 0.20 

1.2 

0.85 

1.1 

2.7 

m-2 

m lo 

= 0.42 

1.2 

0.85 

1.1 

2.7 

m-3 

TTlbr 

= 1.25 

1.2 

0.85 

1.1 

2.7 


Note. — Canonical and Salpeter model ratios, and stellar ratios, included 
for comparison purposes; e Ia = 0.13 and star/gas=0.25 variations of model 
sl-3 also included. 


Table 4. Characteristics of Models with Zp e ,iCM = 0.5 
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iq Ia 
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a SN 
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salpeter 



0.27 

0.11 

0.5 

0.5 

0.25 

0.17 

0.0014 

0.0065 

0.50 

2.0 

0.32 

0.86 

canonical 



0.41 

0.17 

0.5 

0.5 

0.25 

0.17 

0.0022 

0.011 

0.39 

2.0 

0.27 

0.83 

sl-0 

a = 

= 1.83 

0.56 

0.24 

0.96 

1.7 

0.60 

0.22 

0.0020 

0.019 

0.16 

1.0 

0.11 

0.44 

sl-2 

<Uo 

= 1.06 

0.54 

0.22 

0.5 

0.97 

0.42 

0.21 
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0.014 

0.24 

2.0 

0.23 

0.65 
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0.54 

0.23 

0.78 

1.5 

0.54 

0.22 
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0.13 

0.50 

{e Ia = 0.13) 

(Xhi 
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0.45 

0.18 

0.33 

0.75 

0.31 

0.17 

0.0039 

0.013 

0.32 
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0.40 

0.76 

(star/gas=0.25) 
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0.38 

0.16 

0.46 

0.33 

0.43 

0.32 

0.0021 

0.0096 

0.46 

2.2 

0.32 

0.72 
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0.24 

0.93 
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0.24 
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0.65 
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0.54 

0.22 

0.5 

0.98 

0.42 

0.21 
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0.014 

0.23 
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0.23 

0.64 

m-3 
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0.53 

0.22 

0.5 

0.98 

0.42 

0.21 
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0.014 

0.23 
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0.23 
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Note. — Canonical and Salpeter models included for comparison purposes; e Ia = 0.13 and star/gas=0.25 variations of model sl-3 also 
included. 
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mass of star formed is implied by the ICM Fe abundances for otherwise standard assumptions. 
We construct a scenario along these lines by considering an increase in e Ia - the fraction 
of 3 — 8 M 0 that explode as SNIa - from 0.076 to 0.13 (/*/ ficM and stellar abundances 
constrained to match the standard values; Section 4.1). Indeed we find Z Fe j C M = 0.5 
for an IMF that otherwise (i.e., in terms of IMF, star formation efficiency, SNcc lock-up 
fraction and rate, and ICM primordial fraction) approximates the standard model (Tables 
3 and 4). Naturally, the increase in efficiency of formation of SNIa progenitors results in a 
large asymmetry between the stellar and ICM abundance patterns (a SN = 3.1), as reflected 
in the more nearly solar [ a/Fe\icM ratio and consistent with ICM abundance patterns 
(Rjcm = 0-4; Bulbul et al. (2012)) - see Tables 3 and 4, and Figures 9 and 10 that display 
results for such a variation of model si- 3. If this is the correct explanation for the observed 
ICM abundances, one must seek an astrophysical explanation for boosting e Ia in rapidly 
star-forming systems. 


4.4.3. Star-to-ICM Ratio 

As briefly discussed in Section 2.2, the present-day star-to-ICM ratio may exceed our 
standard value of 11%, e.g. due to an unaccounted-for ICL fraction or underestimate of the 
stellar mass-to-light ratio. If we increase /*/ ficM from 0.11 to 0.25, a generally satisfactory 
resolution of the cluster elemental abundance paradox is achieved 9 - a level of Fe enrichment 
and abundance pattern ( Rjcm = 0.32) consistent with observations is attained for an IMF 
and other parameters in line with expected values - with the notable exception of the increase 
in star efficiency to £ s f = 32% - see Tables 3 and 4. The results of this variation on model 
si- 3 are also displayed Figures 9 and 10. It is worth pointing out at this juncture that we 
defined £ s f as the fraction of cluster baryons that form stars; in our models the fraction of 
galactic baryons that form stars (where “galaxies” are defined as locations where stars form 
and eject mass into the ICM) is £ s f/£ ga i = (I + ^gw) - 1 ~ 75% for this /*/ Jicm = 0.25 model, 
but also 2/3 for the standard model. 


5. Summary and Discussion 

Star formation, with a canonical IMF and standard efficiency in producing SNIa, that 
builds up a stellar population comprising ~ 10% of the current overall cluster baryon content 
falls short by a factor of > 2 of explaining a typical rich cluster half-solar ICM Fe abundance 


Stellar abundances are unchanged. 
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(Sections 2. 3-2. 4). This is the case even if predicted ICM abundances are enhanced by 
increasing the efficiency at which metals are ejected from galaxies (and where, as a result, 
the overall abundance in stars is significantly below solar), unless the conversion efficiency 
of cluster baryons into stars is also increased well above 10% (Section 3). 

Section 4 (and Appendix A) constructs and utilizes a phenomenological model for the 
evolution of an old, simple stellar population to quantify the changes in the IMF shape 
(high and low mass slopes, break mass) from its standard form required to bring the ICM 
metallicity and cluster stars into concordance in the sense that they be consistent with the 
same parent star formation history. The necessary departure may be either in the “bottom- 
light” or “top-heavy” sense, with the former tentatively preferred based on better agreement 
with observed ICM abundance patterns and on a higher primordial ICM fraction. It is 
further demonstrated that if a standard IMF is to be preserved, a boost in the efficiency 
of forming stars from gas well beyond that consistent with a gas-to-star ratio of 10, and/or 
of producing SNIa progenitor systems, is required. These calculations are conservative in 
the sense of maximizing the enrichment of the stellar population through the choice of SN 
parameters (e.g., SNIa Fe yields, the upper mass limit for SNcc). 

Stars born in cluster potential wells (or those of their progenitors) must be responsible 
for the high level of enrichment measured in the ICM; however, there is increasing tension 
between this truism and the parsimonious assumption that the stars in the generally old 
populations studied optically emerged from the same formation sites during the same epochs. 
Quantifying this tension, and bolstering the case against the universality of star formation 
are the two primary implications of this study. In the remainder of this section we elaborate 
on these themes. 


5.1. The ICM-Enriching Stellar Population as Distinct from that Observed in 

Elliptical Galaxies 

In some cases the departure from the canonical IMF is modest - a shift of a few tenths 
in the slope over some mass range or an increase in the mass at which the slope steepens from 
0.5 to 1.25 M 0 . However, optical determinations of the IMF from kinematic and population 
studies in elliptical galaxies are trending in the opposite direction (Treu et al. 2010; Auger 
et al. 2010; van Dokkum & Conroy 2010, 2011, 2012a, b; Thomas et al. 2011; Dutton et 
al. 2012a, b; Smith et al. 2012; Cappcllari et al. 2012b; Tortora et al. 2012b; Spinicllo et 
al. 2011, 2012; Ferreras et al. 2012; Sonnenfeld et al. 2012; Goudfrooij & Kruijssen 2013). 
The kinematic evidence that indicates a larger mass-to-light ratio in massive ellipticals than 
expected based on a standard IMF is consistent with either an IMF that is top-heavy and so 



produces more stellar remnants, or one that is bottom heavy and produces more unevolved 
low mass stars. However population synthesis modeling of elliptical galaxy spectra favor the 
latter, and a consensus appears to be emerging for an IMF that is bottom-heavy in elliptical 
galaxies with central velocity dispersions > 150 km s -1 where most of the present-day stellar 
mass in galaxy clusters reside, being at least as steep as a Salpeter IMF if characterized by 
a single slope (and steeper still at the highest galaxy masses). 

The chasm between the amount of metals expected to be produced from a stellar pop- 
ulation with such a steep IMF, and the observed level of cluster enrichment is illustrated in 
the plots of Fe abundance in the ICM versus cluster star-to-gas ratio for three distinct IMFs 
in Figure 11. Standard yield sets and values of the parameters daw (0.5), ./)(“ (0.5), and e Ia 
(0.076) are assumed (Section 4.1). 10 Also plotted are the overall averaged Fe abundances for 
the total cluster baryons; these are independent of the detailed galaxy evolution parameters 
5qw and ff, 0 '- As might easily be inferred from previous considerations, it is clear that the 
model with Salpeter IMF requires an excessively large gas-to-star ratio, and that an IMF 
as steep as 0 ~ m -3 ' 05 unequivocally falls short by more than an order of magnitude of 
producing the required amount of metals. 

The most straightforward explanation for reconciling the steep IMFs based on optical 
spectroscopic studies of elliptical galaxies and the relatively flat IMFs needed to produce the 
cluster metals is to reject the conventional wisdom that the stellar populations in ellipticals 
that dominate the cluster stellar mass are primarily responsible for ICM enrichment. Such a 
decoupling begs the question of the origin and present-day whereabouts of the ICM-enriching 
stars and motivates consideration of scenarios with pre-enrichment (that would presumably 
be accompanied by pre-heating) in protocluster environments by a currently inconspicuous 
stellar population. 

However, an important caveat with respect to the optical spectroscopic studies is their 
general confinement to regions well inside the half-light radius as well as systems at low red- 
shift. Given the emerging paradigm of the multi-stage, inside-out formation/assembly/growth 
of ellipticals by multiple mechanisms (star-forming major mergers, “dry” minor mergers, and 
cold and hot gas accretion; e.g., Conselice et al. 2012, Patel et al. 2012), spatial gradients 
and temporal evolution in properties of elliptical galaxy stellar populations such as the IMF 
is to be expected (La Barbera et al. 2012). With the current dearth of global constraints, 
as well as degeneracies between the inferred dark matter content and the IMF (Wegner et 
al. 2012; Tortora et al. 2012a) and possible systematic errors resulting from nonsolar abun- 
dance ratios (Ferreras et al. 2012), an IMF in cluster ellipticals that is flatter than currently 


10 That is, the solutions are no longer constrained to match the standard stellar abundances. 
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inferred in the core when integrated over space and time is plausible (Dave 2008; Worthey, 
Ingermann, & Serven 2011; Narayanan & Dave 2012b). 

Recent arguments for a reconsideration of bimodal star formation - hints for an IMF 
that is top-heavy in older, but bottom- heavy in younger, star clusters (Zaritsky et al. 2012), 
or top-heavy in denser environments (Marks et ah 2012) - support this. The evolution in 
the cluster star formation environment plausibly leads to an IMF in elliptical galaxies (and 
their progenitors) that transition from one initially weighted towards the mass range that 
includes SNcc and (prompt) SNIa progenitors to one especially conducive to low mass star 
formation at later times (Narayanan & Dave 2012a, b). We note that the high metallicities 
seen in the gas in high redshift quasar hosts (Dietrich et ah 2003a, b) also indicate an early 
epoch of rapid star formation that efficiently produces SN progenitors and is accompanied 
by powerful galactic outflows (Di Matteo et ah 2004; Wang et ah 2002). 


5.2. The Enriching Stellar Population as Distinct from that Observed in our 

Galaxy 

The hypothesis that star formation is universal is refuted by analysis of the level and 
pattern of ICM elemental abundances. The star formation characteristics of the stellar 
population responsible for these metals must depart from that studied locally in one or more 
of the following ways: (1) engender a higher fraction of high mass stars, (2) more efficiently 
form stars from gas, (3) more efficiently produce SNIa progenitor systems. In addition, we 
saw in the previous subsection that the IMF if the enriching population is distinct from that 
recently inferred in the central regions of elliptical galaxy. Arguments for (1) were presented 
above. We now examine the feasibility, and implications, of hypotheses (2) and (3). 

The true star-to-gas ratio (and implied star formation efficiency) remains uncertain, 
with stellar masses difficult to estimate given low surface brightness extended light and 
the likelihood of multiple stellar populations that complicate the conversion from measured 
light in some aperture to total stellar mass (Munshi et al. 2012), as well as the uncertain 
ICL contribution. Both bottom-heavy as now being inferred in elliptical cores, and top- 
heavy/bottom-light as required by ICM enrichment, IMFs may result in upward revisions in 
mass-to-light ratios. A global value of /*/ ficM > 0.25, even for rich clusters, does not seem 
to be excluded by observations at this time. The star formation efficiency corresponding 
to /* / f icm > 0.25, excluding the primordial portion of the ICM that does not engage 
in star formation (see Section 4.4.3, above), is {5e 9a ;(l — r*)} -1 , where e ga i (first defined in 
Appendix A) is the fraction of baryons initially in star-forming structures. This quantity has 
an absolute minimum of 0.2, is > 0.3 for any reasonable value of the mass return fraction 
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r* > 1/3, and > 0.5 for r* >1/3 and e ga i < 0.6. These considerations indicate that a 
substantial fraction of protocluster gas was in the form of dense star-forming protogalaxies 
or pre-galactic fragments. There are clearly profound implications for such a high stellar 
fraction and star formation efficiency for evaluating the magnitude - and perhaps even 
the reality - of the “overcooling” problem (McCarthy et al. 2011) and the physics of star 
formation quenching as it pertains to the cluster environment, as well as for the precision in 
using cluster gas fractions - that must be converted to baryon fractions using a correction 
for stellar content - to constrain the cosmological world model (Allen et al. 2008). 

Given the uncertainty in the nature and possible diversity of SNIa progenitors, and the 
difficulties in reproducing observed rates (e.g., Toonen et al. 2012, Quimby et al. 2012, and 
references therein), the feasibility of an efficiency of SNIa progenitor formation in galaxy 
clusters that exceeds the standard e Ia = 0.076 is not easily evaluated, but cannot be sum- 
marily dismissed. Recent work in this area provides hints, on the one hand, of a downward 
revision in the global estimate of ?/“; but, on the other, of a higher value in galaxy clusters 
(Perrett 2012; Maoz et al. 2012; Graur & Maoz 2012; Quimby et al. 2012). Both an IMF 
that produces additional stars in the 3 — 8 M 0 range, and an increase in e Ia , may boost the 
value of i] Ia . 


5.3. Future Directions 

Progress in resolving the cluster elemental abundance paradox will proceed, in parallel, 
along theoretical and observational lines as follows. Since there is data on the spatial dis- 
tribution of stars and gas (Battaglia et al. 2012), and on the evolution of the Fe abundance 
(Baldi et al. 2012), we are extending our modeling to multi-zone and time-dependent treat- 
ments - with particular attention to possible mechanisms of pre-enrichment and predictions 
for cluster SN (and y-ray burst rates; see below) as a function of redshift and - that fur- 
ther constrain enrichment scenarios. We will also extend our investigation to galaxy groups, 
including fossil groups. 

SN surveys are attaining better statistics, particularly at high redshift, and are sharp- 
ening the accuracy of SN rates, delay-time distributions, and environmental dependencies. 
Optical spectroscopic studies are improving both observationally, and in terms of the com- 
plexity of the stellar population models used to interpret them. X-ray studies of ellipti- 
cal galaxy interstellar and circumstellar gas provide additional probes of elliptical galaxy 
evolution (Loewenstein & Davis 2010, 2012). Future improvements in measuring cluster 
abundance patterns beyond Fe, and in abundance and abundance pattern gradients and 
time-variation are crucial. 
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Finally, we note that many of the mechanisms suggested here for explaining the level 
of ICM enrichment - pre-enrichment by massive stars, efficient and rapid conversion of stars 
to gas, an IMF skewed to high masses - would suggest that the protocluster environment is 
a fertile one for producing y-ray bursts (Lloyd-Ronninget ah 2002; Wang & Dai 2011; Elliot 
et ah 2012), a suggestion we are following up on. 


5.4. Concluding Remarks 

The goal of this work is to quantify the requirements for the stellar population responsi- 
ble for injecting metals into the ICM, and evaluate the feasibility that the stars we see today 
originate from the same source. One is driven to conclude that there is a profound diver- 
gence between the ICM-enriching population and that in the ensemble of elliptical galaxies 
based on standard assumptions about, and recent optical spectroscopic population studies 
of, the latter. This is inferred from the number of SN progenitors needed for the former and 
that expected in ellipticals based on their integrated light and apparent bottom heavy IMF, 
implying the existence of a distinct “hidden” stellar source of metals that may or not inhabit 
the same space as these galactic stars at the same time. While the modeling here is basic, 
the conclusion depend mostly on simple accounting of metals and unlikely to be altered in 
more sophisticated treatments. And although the rich galaxy clusters we consider repre- 
sent an extreme environment, there are broader implications for ellipticals, since mass is a 
much stronger determinant of their formation than environment (Grutzbauch et al. 2011a, b). 
However in it is in the ICM where these phenomena are embedded and remain accessible, 
given the dominance by elliptical galaxies of cluster light, and the closed-box nature of these 
deepest of potential wells. 

We present compelling evidence for a diversity of star formation in terms of some combi- 
nation of efficiency, IMF, and ability to produce SNIa progenitors. Implications to be further 
explored include possible impacts on using cluster baryon fractions to constrain cosmology, 
converting stellar light to mass, and treating star formation and pre-heating/feedback - and 
evaluating overcooling - in semi-analytic models of galaxy formation. 

Occam’s razor is violated in rich galaxy clusters - although metals are made in stars and 
most of the stars we observed are in elliptical galaxies, this stellar population as currently 
understood is evidently not responsible for producing the metals in the ICM. Moreover, the 
nature of the star formation that did produce these metals is clearly very different from that 
we are most familar with, as well as that recently inferred in elliptical galaxies. 


The author wishes to acknowledge useful (and enjoyable) discussions with Richard 
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Mushotzky and Esra Bulbul. This paper is dedicated to my late father, Jerry Loewenstein, 
for reasons that cannot be articulated. 


A. A Simple Model for the Composite Chemical Evolution of Cluster Galaxies 


We approximate the stellar population responsible for enriching the ICM as originating 
from a single, brief, and early star formation episode. As discussed in Section 4, we adopt 
a three-part, continuous, monotonically decreasing, piece- wise power-law form for the initial 
mass function of forming stars (IMF) c f>(m ) = dN/dm following Kroupa (2001); Kroupa 
et al. (2012) extending from mi 0 to rri} n and normalized so that J mfu drn4>(m) = 1. Thus 
4> = kf(x ) for x = m/mi 0 , where 

i 

f(x) = x~ ai Cj, Xi < X < x i+ i, i = 1,3; (Al) 

l=i 


X\ = Ci = 1 , Cj = Cj-iX ^ 3 “ J_1 , X 4 = rrihi/mi 0 , and k is determined from the normalization 
condition. Thus, in addition to the lower and upper mass limits, three slopes (ai,o; 2 , and 
c * 3 ) and two break-masses (m 2 = mi 0 x 2 and m 3 = mi 0 x 3 ) must be specified. 

For the canonical IMF (Kroupa et al. 2012) that we adopt as default, ( mi 0 , m 2 , m 3 , rrihi) = 
(0.07,0.5,1.0,150), where all masses are in M 0 and (an, « 2 ,« 3 ) = (l-3,2.3,2.3). We consider 
IMFs where we vary either m ; 0 (over [0.01,0.5]), m 3 ([0.5,8]); or ai, a 2 , or a 3 (all over 
[0.3, 00 )). 


The mass return fraction for intermediate mass stars is given by 

f mcc dm(fi(rn)Am(rn) 

im _ J >n to r ' ; v ' 

f™ hl dm(j)(m)m 


(A2) 


where rn to = 0.9 M 0 is the main sequence turn-off mass, m cc = 8 M 0 is the lower mass limit 
for SNcc progenitors that we use to delineate “intermediate” and “high” mass stars, and 


A m(m) = m- m rern , wd {m ), 


(A3) 


where 

ITtrem,wd = 0.394 + 0.199m (A4) 

is the white dwarf remnant mass (Kalirai et al. 2008). Similarly, for high mass stars 


hi f 'n UF dm 4 >(m)Am(rn ) 

drri(j)(m)m 


(A5) 



where A m(m) is now based on an averaged remnant mass (see below), 


A m(m) = m- m re m,SNcc{m). 

The specihc number of SNcc explosions per star formed is 

f mup dm(j)(m ) 

cc _ Jm cc 1 

J '™ hi dm(j)(m)m 


(A6) 


(A7) 


A simple chemical evolution model for the composite stellar population in clusters is 
constructed and adopted. We use this to calculate the total mass return from high mass stars 
as the metallicity of the progenitor stellar population is built up, and to connect the ICM 
enrichment parameters with the astrophysics of the formation of cluster galaxies. The model 
is appropriate for stellar populations where conversion of gas to stars is relatively rapid and 
efficient, and so may be applied to cluster galaxies where ellipticals dominate the stellar 
mass and star formation is accelerated in general due to the high primordial overdensity. 
As such, galaxy evolution is divided into two epochs: active and passive, and three phases 
(Loewenstein 2006): star-forming gas (“ISM”), stars, and non-star-forming gas (“ICM”). 
Note that any hot halo gas - relatively insignificant in mass compared to the “true” ICM for 
rich clusters - is subsumed under the ICM category. In the active phase, all star formation 
and SNcc explosions occur and all the initial (ISM) mass in galaxies is consumed by star 
formation or ejected by galactic winds. In the passive phase, stellar mass return and delayed 
SNIa continue to enrich the ICM. 


The cluster as a whole is treated as a closed box, with mass and metal exchange among 
the phases and metal production by the stellar component. We model the active phase 
essentially following the prescription of Qian & Wasserburg (2012) for the case of no infall. 
Mass return is neglected in their approach, and we correct for this in the passive phase. 
Respectively, the evolution equations for the mass in stars, ISM, and ICM are as follows: 


dM* 

dt 


Msf , 


(AS) 


ddVIisM 

dt 


—Msf — M GW , 


dMjcM 

dt 


Mew- 


(A9) 

(A10) 


It is assumed that the rate of outflow is proportional to the rate of star formation that is, in 
turn, proportional to the ISM mass: Msf = A sfMjsm , Mg w = &gwM sf = A sfSgwMism, 
so that the solutions to equations (A8)-(A10) are 


Mjsm — M GALO e xt , 


(All) 
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M* — Mgalo(1 + &gw) 1 (l — e A< ), (A12) 

and 

Micm — Mjcmo + Mgalo^gvk(1 + &gw) : (1 — e A< ), (A13) 

where A = XsfY + Saw), and initial conditions correspond to masses of Mqalo in star 
forming gas (presumably, mostly in galaxies) and Micmo in the ICM, and M* = 0. 

The evolution equations of the corresponding “ith” element mass fractions are as follows: 


dfiJSM _ M^sn 

dt Mjsm ’ 


(A14) 


and 


dfj ,* _ M S f 
dt M * 

dfi,icM _ Mqw 


( fi,ISM ~ /?,*), 

(fi ,ISM — f i,ICM)- 


(A15) 


(A16) 


dt Mjcm 

M it sN (the only source term in the set of equations) is the nucleosynthetic production of the 
ith element, 11 

Mi,sn = 8i,SNMsF, (A17) 

Si,SN V cc (y- c ) + fl a ri Ia yl a , (A18) 


where, as previously defined (see Section 2), y- a and (y™) are the yields per SNIa and SNcc, 
rf c and r] Ia the numbers of SNcc and SNIa explosions per star formed; and, / la is the SNIa 
fraction considered “prompt” in the sense that they occur during the star formation epoch 
(not necessarily part of a distinct prompt SNIa mode). 


Analytic solutions for the metal mass fractions are as follows: 


fijSM = A,svA*f (A19) 

fi,* — $i,SN{ 1 + d G w)~ l t (A20) 

1 — e M 

and 

fijCM = di,SNdcw(l + dcwY 2 , , [l — e~ xt (l + At)] , (A21) 

tviicMyt) 

assuming negligible pre-enrichment of any phase. 

Masses and metallicities at the end of the active phase are assigned according to the 
t — » oo limit of equations (A10)-(A12),(A19)-(A21), following the presumption that most 


11 Only elements primarily synthesized by massive stars and SNIa are considered here. 
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star formation in clusters occurs on a timescale much shorter than the current cluster age, 
and then adjusted for the ensuing passive injection of stellar mass return and “delayed” 
SNIa. The final stellar mass and abundances (including remnants) are, therefore, given by 

m* = m galo ( i - 0(1 + s GW )-\ A* = (i + 5 GW )- l A cc (yT) + % a v Ia y( a ), (A22) 

and the final ICM mass and abundances (elemental mass fractions) by 

Micm = Micmo + MgaloA + $gw) 1 {&gw + r *) (A23) 

and 

fi,ICM = /o(l T &GW ) : [1 + $GW + fo($GW + O] 1 {AgW + ^ *)V° C (vD 

+ [(1 + Saw) - (1 - r*)f v a W a yl a } (A24) 

where the total mass return fraction is r* = r** + r( m . For the overall baryon metallicity 

fiMr = (A25) 

where the star formation efficiency defined in equation (1) is related to e ga i, the “galaxy 
formation efficiency” (the fraction of baryons initially in star-forming structures) according 
to 

£ sf = £ gal(l + 8 gw) \ (A26) 

where e ga i = / 0 /(l + fo) arid f 0 = M galo /M ICM o- 

This formalism enables us to interpret the supernova lock-up parameters introduced in 
Section 2.4 in the context of the chemical evolution of clusters galaxies and place them on a 
firmer physical footing. Naturally the supernova asymmetry parameter a SN = (fp a ) i.e., 
it is the inverse of the fraction of SNIa that explode during the epoch when cluster stars 
form. Our Sections 3 and 4 default aA A = 2 is consistent with estimates of the prompt 
SNIa fraction (Maoz et al. 2011, and references therein; Grauer & Maoz 2012). The lock-up 
fraction may be expressed as /3 CC = (1 — 'r*)(l + 5 gw) _ 1 , and naturally depends both on how 
efficiently stars lose mass, and how efficiently star formation induces galactic winds. 


A.l. Mass Return for Massive Stars 

The active phase chemical evolution model is utilized to calculate the total mass return 
from high mass stars - which can be substantial for top-heavy IMFs. This is motivated by 
the profound impact of metallicity on mass loss in massive stars and, hence, fallback and 
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final remnant mass (Woosley & Heger 2002; Nomoto et al. 2006; Zhang et al. 2008; Fryer 
et al. 2012). The distribution of forming stars as a function of time and mass (Qian & 
Wasserburg 2012) is 


d 2 N Msf 

dmdt dmcj){m)m' 


(A27) 


from which it follows from the expression for the metallicity of star forming gas, equation 
(A18), that 


d 2 N d 2 N (df Fl 


d^tot „—f FB /f Fe n 


0(m), 


(A28) 


dmdfpe dmdt \ dt ) f Fe 0 

where f Fe0 = 5 Fe ,SN/(^ + Saw)- This distribution is used to calculate the mass return from 
massive stars through numerically computing the average remnant mass: 


TFlrem,SNce 


d?N 


Jo 1 dfFe JT;; drn- dmdfFe 


^ j rern,SNcc{^^ / ') f Fe 


d 2 N 


Id d fFe f ZT dm dmdf Fe 


(A29) 


where the remnant mass as a function of mass and metallicity, m rem ^Ncc( m J fFe ), is adapted 
from (Fryer et al. 2012) (delayed explosion scenario) using Fe as a proxy for metallicity. 


A. 2. IMF-dependence of SNIa Rate 


The number of SNIa explosions per star is expected to vary with IMF as follows: 


V 


la 


= 8 


la 


/™//oM>(m) 

lo 

jZZ dm(j)(m)m ’ 


(A30) 


where = 8 M 0 , = 3 M 0 , and e Ia = 0.076 yields the observationally estimated 

fraction of 3-8 M 0 stars that explode as SNIa (Maoz & Mannucci 2012). 
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Fig. 6. — Clockwise from upper left (a)-(d): ICM Fe abundance and Mg/Fe, Ni/Fe, and Si/Fe abundance 
ratios for models with departures, with respect to the standard IMF, in either high- or low-mass slope: 
Aa = a hi — 2.3 for top-heavy/light (sl-1, si- 2) and ai 0 — 1.3 for bottom-light/heavy (sl-3, si- 4 , sl-4), models 
- see text and Table 2 for details. Curves for models si- 3 and sl-4 closely trace each other in the ratio plots. 
The solid and broken horizontal lines, respectively, show the stellar, and standard model ICM, values (Table 
!)• 
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Fig. 7. — Same as figure 6 for departures, with respect to the standard IMF, in lower mass limit or break 
mass - see text and Table 2 for details. 
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Fig. 8. — ICM abundance ratios, with respect to Fe, for Mg - left panel (a), Si - middle panel (b), and 
Ni - right panel (c) - vs. Fe abundance for top-heavy models sl-0, sl-3, sl-f , and sl-5 ; and, for bottom-light 
models si- 2 and m-3 (all bottom-light models considered essentially follow the same curve). 
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Fig. 9. — ICM Fe abundance and Mg/Fe, Ni/Fe, and Si/Fe abundance ratios versus slope above 0.5 M 0 
(assuming slope 1.3 for m < 1 M 0 ) i.e. , model sl-3, for boosted SNIa progenitor formation efficiency ( e Ia ) 
or star-to-ICM ratio (/*///cm)- As for all models, stellar abundances (solid horizontal lines) are fixed at 
their standard values (Table 1). Results for a standard IMF (broken horizontal lines), and for model sl-3 
with standard e Ia and /*/ fiCM, are reproduced. 
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Fig. 10. — ICM abundance ratios with respect to Fe for Mg - left panel (a), Si - middle panel (b), and Ni 
- right panel (c) - vs. Fe abundance for model sl-3 and its boosted e Ia and star-to-ICM ratio counterparts. 
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Fig. 11. — Fe abundance in the ICM (blue curves) and, overall, in cluster baryons (red curves) for the 
following three IMFs: (solid curves) standard (Kroupa et al. 2012), (dotted curves) single slope with (j) ~ 
to - 2 ' 35 (Salpeter 1955), (dashed curves) single slope with <j> ~ to -3 ' 05 . 


